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At very high densities, QCD is in the color-flavor locked phase, which is a color-superconducting 
phase. The diquark condensates break chiral symmetry in the same way as it is broken in vacuum 
QCD and gives rise to an octet of pseudo-Goldstone bosons and a superfluid mode. The lightest 
of these are the charged and neutral kaons. For energies below the superconducting gap, the kaons 
are described by an 0(2) x 0(2)-symmetric effective scalar field theory with chemical potentials. 
We use this effective theory to study Bose-condensation of kaons and their properties as functions 
of the temperature and the chemical potentials. We use the 2-particle irreducible effective action 
formalism in the Hartree approximation. The renormalization of the gap equations and the effective 
potential is studied in detail and we show that the counterterms are independent of temperature 
and chemical potentials. We determine the phase diagram and the medium-dependent quasiparticle 
masses. It is shown that the Goldstone theorem is satisfied to a very good approximation. The 
effects of imposing electric charge neutrality is examined as well. 
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I. INTRODUCTION 

There has been a large effort in recent years to map 
out the phase diagram of QCD as a function of tempera- 
ture and baryon chemical potential [H, [3, d, 0, H, H, 0, [1] ■ 
For example, much work has been done at high baryon 
density and the understanding of this part of the phase 
diagram has improved significantly as compared to one 
or two decades ago. At sufficiently high density and low 
temperature, we know that QCD is in the color-flavor 
locked (CFL) phase 0,0, IQ- This state is a color su- 
perconducting state since the quarks form Cooper pairs 
as electrons in an ordinary superconductor. The attrac- 
tion between the quarks, which renders the Fermi sur- 
face unstable against the formation of Cooper pairs, is 
provided by one-gluon exchange. 

At asymptotically high densities, one can ignore 
the strange-quark mass and quarks of all three col- 
ors and all three flavors participate in a symmetric 
manner in the pairing. The original symmetry group 
SU{3)c X SU{3)l X SU{3)r x U{1)b is broken down to 
SU{3)c+L+R which is a linear combination of the genera- 
tors of the original group. This linear combination locks 
rotation in color space with rotations in flavor space and 
this has given the name to the phase. In the CFL phase 
there is an octet of Goldstone modes which arises from 
the breakdown of chiral symmetry and a singlet arising 
from the breakdown of the baryon- number group U{1)b- 
The latter is a superfluid mode since it is responsible for 
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the superfluidity in the CFL phase. This is analogous 
to the superfluidity encountered in Bose-Einstein con- 
densed phases in condensed-matter systems. Since the 
symmetry-breaking pattern is the same as in vacuum 
QCD, the low-energy properties of the CFL phase can 
be described in terms of an effective chiral Lagrangian 
for the octet of (pseudo)-Goldstone modes and the su- 
perfluid mode [1, M M M 

[l3| . An important differ- 
ence between chiral perturbation theory in the vacuum 
and in the CFL phase is that the latter is at high den- 
sity and the Lagrangian is therefore coupled to chemical 
potentials via the zeroth component of a " gauge field" . 

At asymptotically high densities, all nine modes are 
exactly massless since one can ignore the quark masses. 
At moderate densities, this is no longer the case. The 
quark masses can not be neglected and chiral symmetry 
is explicitly broken. This implies that only the super- 
fluid mode is exactly massless, while the other mesonic 
modes acquire masses. This is relevant for the interior 
of a neutron star. In this case, the quark chemical po- 
tential is of the order of 500 MeV, while the strange 
quark mass is somewhere between the current quark 
mass of approximately 100 MeV and the constituent 
quark mass of approximately 500 MeV [ij]. The mass 
spectrum in the CFL is the opposite of vacuum QCD 
and the lightest massive modes are expected to be the 
charged and neutral kaons K^/K" and /K^. Fur- 
thermore if the chemical potential associated with one 
of the bosons is larger than its vacuum mass, it will Bose 
condense. Various aspects of condensation of kaons in 
the CFL phase have been studied within Nambu-Jona- 
Lasinio models [H, [3 [13, [IE H S [Sll, [H ■ 

Some of the properties of the kaons in the CFL phase 
have been studied using effective scalar field theories. 



For example, it has been shown that the symmetry- 
breaking that accompanies Bose condensation of kaons 
gives rise to unconventional Goldstone bosons if the orig- 
inal symmetry group is SU{2) x U{1) [H, [13, Hi] . This 
is the case when the chemical potential for the neutral 
kaons is the same as that of the charged kaons. The 
Goldstone mode has a quadratic dispersion relation for 
small values of the three-momentum instead of the usual 
linear dispersion relation. This unsual property of the 
Goldstone mode arises from the lack of Lorentz invari- 
ance due to finite density. 

In contrast to hadronic matter in heavy-ion collisions, 
bulk matter in compact stars must (on average) be elec- 
trically neutral and so a neutrality constraint must be 
imposed [H, H^l • Similarly, bulk matter must be color 
neutral and if the system is in a color superconducting 
phase, one sometimes has to impose this constraint ex- 
plicitly. It is automatically satisfied if one uses the QCD 
Lagrangian, but not so if one describes the system using 
NJL-type models. If there are dynamical gauge fields 
present, the zeroth components of A° develop nonzero 
expectation values, (Ag) ^ 0, that effectively act as 
chemical potentials. In the NJL models there are no 
gauge fields present and the SU{Nc) color symmetry is 
global m, m, m, HOl . One must therefore impose these 
constraints explicitly. 

One of the most popular approaches to the study 
of systems at finite temperature and density is the 
2-particle irreducible (2PI) effective action formalism 
developed by Cornwall, Jackiw and Tomboulis [3l| 
in the context of relativistic field theories. See e.g. 
Refs. Hi, m H m, m, [13, m, m . for various applica- 
tions. The 2PI functional depends on background fields 
(j)o and the exact propagator D. In practical calculations 
one must truncate the exact generating 2PI functional 
according to some approximation scheme. The 
expansion and the loop expansion are two such system- 
atic schemes. These approximation schemes are nonper- 
turbative in the sense that they sum up loop diagrams 
from all orders in the loop expansion. The renormal- 
ization of such approximations is therefore a nontrivial 
issues as standard theorems from perturbation theory do 
not apply. There has been significant progress regard- 
ing the renormalization of these approximations start- 
ing with the papers by van Hees and Knoll [40]. They 
showed that the equation for the two-point function in 
scalar (/)^-theory can be renormalized by introducing a 
finite number of counterterms. A more systematic study 
of the renormalization issues in scalar 0'*-theory was pre- 
sented in Ref. [4l| . They formulated an iterative renor- 
malization procedure that determines the counterterms 
needed to eliminate the (sub)divergences. These coun- 
terterms were shown to be independent of temperature. 
In Ref. [12] the authors have developed a "direct" renor- 
malization procedure which is equivalent to the iterative 
procedure in [4l|. Later it has been shown ^] how to 



fix all the counterterms needed to calculate the proper 
vertices which are encoded in the effective action. These 
counterterms are local and they are independent of tem- 
perature and chemical potential. 

The 2PI effective action formalism for scalar fields has 
recently been used to study the thermodynamics of pi- 
ous and kaons and their condensation. In Ref. ^39], the 
quasi-particle masses and the phase diagram is stud- 
ied to leading order in the expansion of 0{N)- 
symmetric scalar field theories. In Ref. the au- 
thors applied the 2PI effective action formalism in the 
Hartree approximation to an effective 0{2) x 0(2)- 
symmetric scalar field and calculated the phase diagram 
and the critical temperature for Bose-condensation of 
kaons. The effects of imposing electric charge neutral- 
ity were also investigated. The scalar theory for the 
kaons were derived from the effective chiral Lagrangian, 
where the parameters depend on the baryon chemical 
potential. Renormalization issues were not addressed. 
In the present paper, we reconsider the problem of kaon 
condensation from a somewhat different angle, and con- 
sider in some detail the renormalization of the theory. 
We also show that the violation of Goldstone's theorem 
is negligible. 

The paper is organized as follows. In Sec. II, we 
discuss the 2PI Hartree approximation for a 0{N)- 
symmetric Bose gas. In Sec. Ill, we briefiy discuss 
0(2) X 0(2)-symmetric models, which are relevant for 
kaon condensation in the CFL phase of dense quark mat- 
ter. We determine the quasi-particle masses as well as 
the phase diagram. Finally, we study the effects of im- 
posing electric charge neutrality. In Sec. IV we summa- 
rize and conclude. In appendix A, we discuss the renor- 
maliztion of the gap equations and effective potential in 
detail. 



II. 0(2Ar)-SYMMETRIC BOSE GAS 

The Euclidean Lagrangian for a Bose gas with N 
species of massive charged scalars is 

^ - (9^'i>l)(5^<f.) + "^'*l*. + ^(<&l<i>,)' , (1) 

where i — 1,2,.. .,N and — (02i-i + 'i4>2i)/V2 is 
a complex field. The theory described by Eq. ([1]) has 
(2iV — 1)N conserved charges which equals the number 
of generators of the group 0{2N). 

A gas of TV species of bosons can be characterized by 
the expectation values of the different conserved charges 
in addition to the temperature. For each conserved 
charge Qi, one may introduce a nonzero chemical po- 
tential ^i. However, it is possible to speciiy the expec- 
tation values of different charges only if they commute. 
The maximum number of commuting charges is N (43 | 
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and these can be chosen as 



where the the current densities j'j^ are 
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2i 
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(3) 



The incorporation of a conserved charge Qi is done by 
making the substitution 



(4) 
(5) 



in the Lagrangian ([T|) . Note that the chemical potential 
acts as the zeroth component of a gauge field. 

From the path-integral representation of the thermo- 
dynamic potential Q 



-fivn _ 



-f^'drfd^xC 



(6) 



the expression for the charge density can be written as 
ft = -1^- (7) 



If we introduce k chemical potentials, the full symmetry 
group is broken down to [0(2)]'= x 0(27V-2fc). Ifm^ < 0, 
the 0{2N) symmetry is spontaneously broken down to 
0{2N — 1). Even if > 0, the symmetry may be 
broken if one of the chemical potentials, ^i, is larger than 
a critical chemical potential /Xc = m. In that case, the 
(O|02i-i|O) 7^ and the corresponding 0(2) symmetry 
is broken. In the following, we consider the simplest 
example of a single chemical potential /i — fii for the 
complex field $i. We first introduce a nonzero vacuum 
expectation value 4>o for the field $i in order to allow 
for a charged condensate. Using the 0(2)-symmetry, we 
can always choose (/)o real and so we can write 
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(8) 



where and (f>2 are quantum fluctuating fields. 



The inverse tree-level propagator then reads 
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where w„ = 2nnT are the Matsubara frequencies and ^^2 _ ^2 ^ ^ ^2 ^-|^2) 

the tree-level masses are 2N 



Q 3A , o 

H 

2N 



: (10) The tree- level dispersion relation is found by analytic 

continuation to Minkowski space, ujn — + iuj, and then 
-I- i^'fo J (11) solving the equation DetI?o('^jP) — 0- This yields 
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W3(P) 



+ — {ml + 7712 ) + ± \j 4/i2 



p2 H — {m\ + 7712) 
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(13) 
(14) 



The classical effective potential is 
V 



The minimum of the classical potential V is given by 



(m^ - M^) + A/2iV</«i^ = and at the minimum, we have 
nil — 3/^^ — 2to^ and — /i^. The dispersion relation 
'^1,2 (p) then reduces to 

(16) 

From this equation, we note that uj2{p) is a massless 
mode. Expanding around zero momentum p, we find 



U}2(p) 



(17) 



The mode is linear in the momentum for small momenta 
and is thus a conventional Goldstone mode. This is in 
agreement with the fact that the 0(2)-symmetry has 
been broken due to the condensate (po . 



A. Effective Action and Gap Equations 



The 2PI effective action can be written as 



n\4,o,D] 
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-TrlnD" 



(18) 



where D is the exact propagator and ^[D] is the sum of 
all two-particle irreducible vacuum diagrams. The trace 
is over field indices as well as over space-time. In the 
Hartree approximation one includes the double bubble 
diagrams shown in Fig. [T] If we denote by Dij the com- 
ponents of the propagator D, the functional $[-D] can 
be written as 



Generally, in the vacuum, the terms in the 2PI effective 
action p8)) can be classified according to which order 
in the 1/iV-expansion they contribute and can be ex- 
pressed in terms of 0{N) invariants such as Tr(£'") and 
Tr((/)gi:)") The first and second term in Eq. ^ 

are proportional to Tr[0g] and [Tr(0Q)]^, respectively. 
Each trace gives a factor of N and so they both scale 
as N. The third term is the trace of the full propa- 
gator and therefore scales as N. The third term can 
be decomposed into a sum of the terms Tr[D], Tr[(/)gD], 
Tr[(?!)2]Tr[L»], and Tr[02£)]. Finally, let us consider 
Performing the sums involving Fijki, ^[D] can be ex- 
pressed in terms of the invariants Tr{D^) and [Tr(Z?)]^: 



^ ^ 87V 



[Tr (D)]^ + 2Tr (D^) 



(22) 



The first term involves two traces each giving a factor 
of N. There is a factor of 1/iV coming from the vertex 
and so this term goes like N. Similarly, the second term 
involves one trace and one vertex and it contributes at 
order one. We therefore conclude that the Hartree ap- 
proximation is not a systematic approximation. It is not 
systematic in powers of as we have seen and it is 
not systematic in number of loops since we have not in- 
cluded the setting-sun diagrams that arise in the broken 
phase. 





(19) 



where Fi^^i = (<^ij4i + <5jfe<5jz + 5u6jk) is the sum of the 



three rank-four invariants of 0{N). 
above is defined by 



The sum-integral 



47r 



T 



qo=27TnT 



(2^ 



(20) 



where Q — {qo,q), d = 3 — 2e and A is the renormal- 
ization scale associated with dimensional regularization. 
The sum is over Matsubara frequencies. The integral 
over three-momentum q is calculated with dimensional 
regularization. We also introduce the compact notation 
for these integrals 



(21) 



FIG. 1: Vacuum diagram contributing to the effective po- 
tential. in the Hartree approximation. Left diagrams are of 
order N, while right diagrams are of order one. 



At the stationary points, the 2PI effective action sat 
isfies the gap equations 



SD 



= 
= 



The gap equation 



can be rewritten as 
<5$[D] 



(23) 
(24) 

(25) 



Using the fact that D^^ - D^^ ^ n(P), where n(P) is 
the self-energy, we obtain 

m[D] 



U{P) = 2- 



5D 



(26) 
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Note that the self-energy in Eq. is a matrix. The 
self-energy is obtained by cutting a line in the vacuum 
graphs. A generic loop diagram that contributes in the 
Hartree approximation is shown in Fig. [21 




FIG. 2: Feynman diagram contributing to the self-energy 
n(P) in the Hartree approximation. 



In the Hartree approximation, the self-energies are 
momentum independent. Moreover, it can be shown 
that only the diagonal elements of the self-energy matrix 
Hij are nonzero. If we denote these elements by 11^, the 
dressed inverse propagator is given by Eq. Q with the 
tree level masses mi, m2, ma replaced by the effective 
medium-dependent masses Mi, M2, M3, where 

Mf = ml + n, . (27) 

The full propagator then reads 



D{uJn,P) 



-2flUJn 



n ' 1 



(a;2+cD2)(^2+^2) (^2+^2)(^2+^2) 






V 



(28) 



where tii_2,3(p) are obtained from Eqs. ()13|) and ()14p by the replacement nii 
be written out cxplcitly as 

+4(iV - l)jf D22j{ D33 + ^N{N - l)jf Dssjf D. 

JQ JK JQ JK 



Mi. The functional can then 



D 



33 



K 



(29) 



Notice that the terms involving off-diagonal elements of the full propagator D are absent. This is due to the fact 
the diagram vanishes upon summation over the Matsubara frequencies. This follows immediately from Eq. (|28p . 

The gap equations for the dressed masses now follow from Eqs. (|25p . (|26p . (|27p . and These equations 

contain ultraviolet divergences and require renormalization. Renormalization is discussed in Appendix A and the 
diagrammatic interpretation of the iterative renormalization procedure is shown in Fig. [3l 






FIG. 3: Diagrammatic interpretation of the gap equations for rm in the 2PI Hartree approximation. 
After renormalization, we obtain 

Ml -ml = ^ [3( Jf + Jl) + ( J| + Jl) + {2N - 2){J^ + jf)] , (30) 
Mi - "^2 = M + + 3(^2 + Jl) + {2N - 2)iJ^ + 4)] , (31) 
Mi - "^3 = [("^f + + + Jl) + 2N{J^ + J^)] , (32) 
where the integrals and jj are defined in Appendix A. 
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We next consider the gap equation ([M)) . The expression for it follows from and After renormalization, 
we obtain 



= 



- + + ^ [3(Jr + ) + {J!i + J2^) + (2Af - 2)(J^ + Jl')] 



(33) 



r 



Using the gap equation pO|) to eliminate the integrals 
and J^, we can rewrite Eq. (j33p in a very simple way 



= 



Mt 



2 i,2 



(34) 



associated with the breaking of the 0(2)-symmetry due 
to the condensate 0o- Writing M| — ^'^ + S, where 
^ ~ W [("^2 + "^J) ~ ('^i + Ji)] ' ^ calculation analogous 
to the one leading to Eq. ([T7)) . shows that the mass gap 
of the Goldstone mode is given by 



The difference between this equation and the corre- 
sponding equation at the tree-level, is the replacements 
rrii — > Mf. Similarly, combining the gap equations pip 
and (|33p. we obtain 

= + ^ [iJ^ + Jl) ( + JD] ■ (35) 

Comparing this equation with the tree-level result — 
fi^, Eq. (j35p . we see that the Goldstone theorem is not 
respected. It is well known that the 2PI Hartree ap- 
proximation violates Goldstone's theorem and in the 
present case this means that there is no massless mode 



Ml 



(36) 



In the 2PI 1 /A^-expansion, the Goldstone theorem is sat- 
isfied order by order. In the large- A' limit, it is clear from 
Eq. (|35p that the Goldstone theorem is respected since 
the last term vanishes. 



The effective potential ^ follows from Eqs. ([9]), ^ 
and (j29p It can be renormalized in the same manner as 
the gap equations and details can be found in Appendix 
A. After renormalization, we obtain 



n ^ -{n?^ ^?)ct.l + —c^t + -{Jl + J^) -{J^ + Jl) + {N- \){J^ + Jl 

-\{Ml - niDiJ^ + Jl) - i(Af| - m2)(J| + jj) - (A^ - l){Mi - ml){J^ + 4) 



where J^ and J^ are defined in Appendix A. 

III. KAONS IN THE CFL PHASE 
A. Effective theories 



ventional ipip. Nevertheless, one finds the same quantum 
numbers. Another difference is that Lorentz invariance 
is broken due to the presence of the chemical potentials 
and so the Lagrangian in invariant only under rotations 
in three dimensions. The chiral effective Lagrangian is 
given by 



In the CFL phase of dense QCD, the original symme- 
try SU{3)cX SU{3) L X SU{3)r xU {1)318 broken down to 
SU{3)c+L+R- The diquark condensate {il'tp) breaks chi- 
ral symmetry in exactly the same manner as in vacuum 
QCD and so the effective Lagrangian for the Goldstone 
modes have the same structure as in chiral perturbation 
theory. Notice, however, that the mesons are composed 
of four quark fields of the form iptpipip instead of the con- 



£ = -/2Tr[(5oS + z[AS])(aoS-z[A,S]t) 

-vl{^.m^^^^)] 

+ ^af^detMTr[M-\j: + Y.^)] +■■■, (38) 
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where w^, and a are constants. The meson field S is 
given by 



(39) 



where A° are the Geh-Mann matrices and 0° describe 
the octet of Goldstone bosons. The matrix A acts as 
the zeroth component of a gauge field and is given by 



A 



M2 

2fj. ' 



(40) 



where fig is the chemical potential for electric 
charge, /i is the baryon chemical potential, Q = 
diag(2/3, — 1/3, — 1/3), and M — diag(m„, to^, nis). 

At asymptotically high densities, one can use pertur- 
bative QCD calculations to determine the parameters 

, , and a by matching [l^, [l^, [l^ ■ 



f = 



12-81n2 
187r2 ' 



(41) 



1 

3 ' 
3A^ 



(42) 
(43) 



where A is the superconducting gap. Note that = 
l/\/3 is the standard result for the speed of sound in 
a dense medium. In the vacuum, Lorentz invariance 
enforces the value = 1. It is important to stress 
that the values of these parameters given by Eqs. (|4T|) - 
(|43p are valid only at asympotically high densites. For 
moderate densities, which are relevant for compact stars, 
one does not know the values of these parameters as they 
cannot be determined by matching. 

After expanding to fourth order in the meson fields, 
Alford, Braby, and Schmitt obtain the following effective 
Lagrangian for the kaons: 



C 



(do 
1 



Ml 



Po-Aii)*i 

/32 



(a,$I)(9,$i) + [(ao + Ai2)$^ 



[(^0 - Ai2)$2] + 2"^?*1 



*2 + ^'^I'^'l , 



(44) 



where the complex doublet are {K^, if°) = ($1, $2) and 
the parameters are given by 



m,l = amd{ms + mu) 
m\ = amu{ms + m^) 



Ml = MQ + 



9 2 
mi — mt 



2fi 



M2 



2fi 



1 r« , A /^Mi - M2 



(45) 
(46) 

(47) 
(48) 
(49) 

(50) 



The dimensionless parameters a and Pi are quartic cou- 
plings of the effective theory. Note in particular that 
they depend on the chemical potentials ^i. The effec- 
tive theory described by the Lagrangian (|44| is invariant 
under the group 0(2) x 0(2). 

There is a technical complication arising from using 
the Lagrangian Eq. (|44|) . The problem is that the effec- 
tive couplings depend on the chemical potentials. This 
implies that the counterterms also depend on the chem- 
ical potentials i.e. parameters that describe a dense 



medium. Renormalizing a theory based on the La- 
grangian (|44p therefore depends on the medium, which 
one may object to. We therefore take a somewhat dif- 
ferent approach by using an effective Lagrangian with 
mass parameters and couplings that are independent of 
the chemical potentials. For simplicity, we set f,r = 1, 
but it is not difficult to scale loop momenta in our equa- 
tions to take into account values of Vj^ that differ from 
the its vacuum value. The kaons are written as a com- 
plex doublet \ {K°,K+) = ($i,$2). The Euclidean 
Lagrangian with an 0(2) x 0(2) symmetry is given by 



L = 



(5o + Mo)$l] [(^0 - Mo)$i] + (a,$l)(a,$i) 
(So + /.+)$tl [(ao - ^+)<i>2] + (a.$t)(a^$2) 



-m^$|$i -I- m\<^\^2 + ^ (*I*i 



(51) 



The chemical potentials /io E^nd M+ associated with the 
two conserved charges for each complex field ^i. They 



1 Note that the identification in Ref. [3 is (K+,K°) = {<I>i, 'I>2). 
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are related to the quark chemical potentials fj.u, fJ^d, and 
Ms by 

/^o = y.d- tJ-s , (52) 

= fj-u- fj-s ■ (53) 

We therefore have /i+ — /ig = /iu — /i^- In weak equilib- 
rium, the processes d+v <-> u + e^ go with the same rate 
in both directions. If we assume that the neutrinos leave 
the system, their chemical potential is Hi, = 0. This im- 
plies that iJ,d = l^-u — ^J'Q^ where fig is the electric charge 
chemical potential. In other words, /ig — /ig. In the 
absence of the chemical potentials and with mo = m+ , 
and Ao = A_|_ = A^f, the Lagrangian ([5T|) has an ex- 
tended SO{i) ~ SU{2)lX SU{2)r symmetry If we add 
chemical potentials such that /io = this symmetry 
is broken down to SU{2) x J7(l), where the chemical 
potential is for the C/(l) charge. A condensate would 
then break this symmetry down to U{\) implying the 



existence of massless modes. Naively, one would per- 
haps expect three Goldstone modes as there are three 
broken generators. However, one of the massless mode 
is quadratic in the momentum p for small p and the 
Nielsen-Chadha theorem implies that such a mode be 
counted twice [25|, ,48, J9J . This is consistent with the 
fact that there are only two massless modes. 



B. Effective potential and gap equations 

In order to allow for a condensate of neutral kaons, 
we introduce an expectation value (/jq for $i and write 

$1 = --^(</-o + 0i+i02) , (54) 

where (f)i and 4>2 are quantum fluctuting fields. 



I 

The inverse tree-level propagator can be written as a block-diagonal 4x4 matrix: 



ml 



Mo 








UJ„ 



— 2pQU!n 

-t- m^ 





Mo 






+ m 

2p+UJn 






-2p+UJn 



p- 



J 



(55) 



where the tree-level masses are 

2 



m 



m 





3Ao 




ml + 


2 




ml + 


2 ^ 


I. 




\h 




m\ + 


2 


00 



(56) 
(57) 
(58) 



Note that in the remainder of this section, the mass 
parametei 
tential is 



parameters ml and m^ are positive. The classical po- 



V=\{ml-^?,)4>l + ^4 



The dispersion relations are 
^1,2 (p) = 



\ 



p^ + -{ml+ml)+lll±JA^ll 



p2 + -{ml+ml) 



1 

4 



-\ — {m\ — m^ 



2 

21 I 



The 2PI effective action is given by 

1 



(59) 



(60) 
(61) 

(62) 



where $[-D] contains all 2PI vacuum diagrams. In the Hartree approximation, we include all double-bubble diagrams 
which can be written in terms of 0(2) x 0(2) invariants. If we denote by Da and Dh the two 2x2 submatrices of 
the propagator D, we can write [i^l 

^D] = ^ [Tr {Da)f + 2Tr {dI)] + ^ [Tr {D,)f + 2Tr + ^(TrZ?,)(Tr7?,) . (63) 
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Writing out explicitly the terms in $[-D], we find 
3 



<^[D] = I Aorf Dnj{ Dn + Aorf ^^22^ D22 + X+jf Ds^jf D33 + X+jf D^^jf dJ 
* L Jk Jq jk Jq Jk Jq Jk Jq J 

4 L Tft: Tq Ta' Jq Jk Jq Jk Jq 

+XhJ{ D22j{ D33 +XhJ{ D22j[ Da 
Jk Jq Jk Jq 



(64) 



Again the terms involving the off-diagonal elements of the full propagator D are absent since they vanish upon 
summation over the Matsubara frequencies. 

The gap equations now follow in the usual manner. They require renormalization, which is briefly discussed in 
Appendix A. After renormalization, the gap equations read 



Ml - ml 



Ml - ml 
M| - ml 



\ [3Ao(Jr + Jf ) + Ao(J2^ + Jj) + 2Xn{Jl + JI)\ , 
i [Ao(Jr + Jl) + 3Ao(J| + Jj) + 2Xh{JI + Jl)\ , 
\ {Xh{J{ + Jl) + Aff (J| + Jj) + 4A+(J3^ + Jl)\ , 







^0 - /^o + ^-^0 + \ [3Ao(Jf + Jl) + Ao(J2^ + Jl) + 2Xh[,JI + 4)] 



(65) 
(66) 
(67) 
(68) 



In the appendix we argue that M4 is equal to M3 and their gap equations are idential. This is only correct as long 
as there is no charged kaon condensate. Substituting Eq. ([65t into Eq. ((68)) . it can be written as 



</)o[M2^(Mg + Ao<^2)] ^0, 
Similarly, substituting Eq. ([55]) into Eq. (1551) . we obtain 

Ml -ill = X^ {{Jl + Jl) - {Jl + Jl)\ 



(69) 



(70) 



The effective potential can be renormalized using the same methods and after renormalization the effective po- 
tential is given by 
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1 



1 



{Jl + Jl) + - {J!^ + Jl) + ( + Jl 



\{Ml - m\){Jl -t- Jl) - \{Ml - m\){Jl + jj) - (M| - ml){Jl + J^) 



"^^0 / TC , tT\2 , 3Ao 



J2 ^ "^2 



Ao 



'A + Jl) {J2 + Jl 



(71) 



r 



C. Phase diagram and quasiparticle masses 

In order to determine the phase diagram and numer- 
ical evaluate the masses of the quasiparticle, we need 
to know the values of the parameters in the effective 
Lagrangian ([5TII . In Ref. [IJ], the authors use a quark 
chemical potential [i ~ 500 MeV and and superconduct- 
ing gap A ~ 30 McV. Extrapolating from asymptoti- 
cally high densities, Eq. (gll), and Eqs. (gHll-lgS]) give 



U - 100 MeV, ~ 20 MeV, mi ~ 5 MeV, and m2 ~ 4 
MeV. In the numerical calculatations that we present be- 
low, we use mo = 4 MeV, m+ ~ 5 MeV, /io = /i+ = 4.5 
MeV, and /tt — 100 MeV unless otherwise stated. The 
renormalization scale is chosen to be the average of the 
two vacuum masses, i.e. A = 4.5 MeV. Using Eqs. ([35])- 
(|50p and identifying 2/3i with Aq, 2/32 with A+, 2a 
with A/f, the numerical values are Aq = 1.25 x 10"'^, 
A+ = 1, .08 X 10-3, and Xr = 1.16 x 10"^. In the plots, 
where we vary the chemical potentials, we use the same 
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couplings throughout as /i-dependcnt couphng constants 
are problematic. 

The neutral kaons K'^ and if are identified with the 
linear combinations of the fields 0i and (j)2 , and therefore 
with uJiijp) and Cj2{p). Similarly, the charged kaons K~ 
and are given by linear combinations of 03 and </)4 
and are identified with Cj^(p) and Oji{p). 

In Fig. m we show the neutral kaon condensate as a 
function of /ig and /i+ for T = 0. For /i+ = 0, i.e. 
along the ^g-axis, there is a second-order phase transi- 
tion to a neutral phase with a kaon condensate at a crit- 
ical chemical potential /io — m2- This is the CFL-if" 
phase. For larger values of /i+ the transition becomes 
first order. The point in the (/ig, /i+)-plane where the 
transition changes order is a critical point and given by 
(4.0, 5.0) MeV. In the part of the phase diagram where 
the transition is first order, the transition is actually to 
a phase with condensate of charged kaons. This con- 
densate is not shown in the figure. This is the CFL-i^+ 
phase. Thus there is a competition between the neutral 
and the charged condensates and nowhere do they exist 
simultaneously. The transitions to the kaon-condensed 
phases are density driven. 




FIG. 4: Neutral kaon condensate as a function of the chem- 
ical potentials /io and ^+ for T = Q. 



In Fig. [5l we show the neutral kaon condensate as a 
function of /io and /i+ for T = 200 MeV. The point in 
the (/xq, /i+)-plane where the transition changes order is 
given by (5.2,6.0) MeV. 

In Fig. [HI we show the mass parameters Mi^2 and 
normalized to /io for /io = //+ = 4.5 MeV as a function 
of T normalized to Tc- The masses M3 and M4 are de- 
generate for all values of T, while Mi and M2 become 
degenerate at the critical temperarture. If the Gold- 
stone theorem is obeyed M2 is exactly equal to in the 
broken phase. We notice that there is a tiny deviation. 

In Fig [71 we show the thermodynamic potential 
^{4>o) ~ ^{0) as a function of the condensate 0o for 
/io = /i+ = 4.5 MeV and three different values of the 
temperature. fl{(f>o) is obtained by solving the gap equa- 




FIG. 5: Neutral kaon condensate as a function of the chem- 
ical potentials fio and /i+ for T — 200 MeV. 




0.0 0.5 1.0 1.5 



FIG. 6: Dressed masses Mi, 2 and M3 normalized to /io as a 
function of T/Tc for /to = /i+ = 4.5 MeV. 



tion for the masses (pS)) . and ((55)) and inserting the 
values into the effectifve potential ([7T|) . The solid line 
is T = 0, the dashed line is T = = 118.5 MeV, and 
the dotted line is T = 200 MeV. The phase transition is 
second order. 

In Fig. [5] we show the masses of (1^2 (g = 0)) and 
{uj4{q — 0)) for ^0 = M-i- = 4.5 MeV and as func- 
tions of T normalized to Tc. We notice that the mass 
of K'^ is not strictly zero, which explicitly shows that 
the Goldstone theorem is not respected by the Hartree 
approximation ^. In Ref. (T^ . the authors make some 
further approximations of the sum-integrals appearing 
in the gap equations. These approximations give rise to 
an exactly gapless mode. As pointed out in their pa- 
per and as can be seen in Fig. [51 this is a very good 



^ Note that the mass gap does vanish at T = T^, since Jj = 
and Jf = Jj, cf. Eq ^Si- 
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40 60 
00 [MeV] 



FIG. 7: n{4>Q) - fl{0) for fio ^ n+ ^ 4.5 MeV and three 
different values of tfie temperature. The sohd line is T = 0, 
the dashed line is T = Tc = 118.5 MeV, and the dotted line 
is T = 200 MeV. 



approximation. 
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FIG. 8: Mass gaps for the K'^ and Tf" modes for fio = 
4.5 MeV and as a function of T normalized to Tc- 



In Fig. [HI we show the difference M| - fj.^, Eq. (|7D)) . 
normalized to for fiQ = 4.5 MeV as a function oiT/Tc- 
This is a another measure of the violation of Goldstone's 
theorem. We see that the violation is tiny. 



D. Effects of charge neutrality 

Bulk matter in compact stars must be overall color 
and electrically neutral, otherwise one pays an enor- 
mous energy penalty. The neutrality constraint applies 
whether or not the gauge charges are broken or not. For 
certain values of /Lto, and T, the thermodynamically 
stable state, i.e. n evaluated at the stationary points, 
has an overall electric charge. For example, in the re- 
gion of the lower corner of Fig. [H there is a charged kaon 
condensate which is positively charged. 

In this subsection, we impose the constraint of overall 
charge neutrality. In some case this constraint can have 



o 
3, 



3. 
I 




0.0000 



FIG. 9: Measure of the violation of Goldstone's theorem, 
(M| - Mo)/^*o as a function of T/Tc. 



dramatic effects on the phase diagram. For example, 
in a recent study by Abuki et al [50|, they showed that 
pion condensation in the presence of finite isospin chem- 
ical potential does not occur for physical masses of the 
pions. In fact, they found a tiny window of pion con- 
densation for pion masses below approximately 10 KeV. 
Thus pion condensation is very sensitive to the explicit 
chiral symmetry breaking. 

For simplicity, we assume that we can describe the 
background of electrons by an ideal Fermi gas. We then 
add to the Lagrangian ([5T|) . the term 



electrons 



(7^5^ + 7Vee - me) Ve , (72) 



where i/jg denotes the electron field, e is the electron 
charge, and me is the mass of the electron. The contri- 
bution to the free energy from the electrons is denoted 
by r^e and reads 



r^e = -2 

-Tlog 



j {^p + Tlog 



I _ ^-P(Ep-l-iQ) 



,-/3(£;p+MQ) 



(73) 



where Ep — y/jp 



In the following, we assume the 



electrons are massless. If we use dimensional regular- 
ization, the first term in Eq. (I73p vanishes since there 
is no mass scale. The temperature-dependent integrals 
can be done analytically and we obtain 



2^Q 



6 



2 ji2 



111 

180^ 



(74) 



We add Eq. to Eq. ([IT]) to obtain the full thermo- 
dynamic potential. The contribution to the electric 
charge density from the kaons is then given by 



71+ = 



(75) 



11 



Using Eqs. and ([TS]), we obtain 



n+ 



(76) 



This expression is free of ultraviolet divergences. After 
summing over Matsubara frequencies and averaging over 
angles, the equation reduces to 



1 

4^ 



1 



1 



dq 



(77) 



The contribution He to the electric charge density from 
the electrons is given by 

Charge neutrality amounts to requiring that 

n++ne . (79) 

In the case where we do not impose the charge neu- 
trality requirement (|79p . we have two chemical poten- 
tials that we can vary freely. The charge neutrality con- 
straint gives a relation between /io and /i+ and only one 
of them is free to vary. For a given value of T and 
e.g. Ho, we must therefore solve simultaneously the gap 
equations ((65)) - ((70| and the equation (|79|) . to find the 
dressed masses Mi, 0o and the chemical potential In 
Fig. [ini we show the neutral kaon condensate for neutral 
matter as a function of T and /iq For T = 0, the onset 
of kaon condensation is at /iq = ruQ as expected. The 
transition is second order for all values of /Lto and again 
the transition is to a symmetric phase. 



the curves are always equal to one. The chemical po- 
tential — /iQ vanishes at T = 0, increases rapidly, and is 
essentially constant for T > 20 MeV. The chemical po- 
tential which equals /xq at T = 0, never reaches its 
critical value /x'^ = m_|_ and this explains why the phase 
transition in Fig. [10] is to a symmetric phase and not to 
a charged kaon condensed phase. In other words, charge 
neutrality does not allow for a charged condensate. 
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FIG. 11: The chemical potentials and — /iq divided by /io 
as functions of temperature in an electrically neutral system. 



Finally we compare the critical temperature with and 
without neutrality. In the case with no neutrality con- 
straint, the critical temperature was 118.5 MeV for 
/io = /i+ = 4.5 MeV. In the neutral case, /i+ is a func- 
tion of /io and so varies with temperature. At T = 0, we 
always have /xg = 0, i. e. /io = M-i-. We therefore choose 
the same value, i.e. /io — 4.5 MeV. The critical tem- 
peraature in this case is = 125 MeV. In other words, 
there is a small increase in the critical temperature as 
one impose charge neutrality. This is in agreement with 
the findings of Ref. [11]. 



IV. SUMMARY AND OUTLOOK 



2m ■ 



\m 



^1 



FIG. 10: The neutral kaon condensate as a function of T and 
/io for a system which is electrically neutral. 



In Fig. 1111 we show /i_(- and — /ig = /io — /i+ normal- 
ized to /xo as functions of temperature T. The sum of 



In the present paper, we have studied the 0{N) and 
0(2) X 0{2) models with chemical potentials using the 
2PI formalism in the Hartree approximation. We have 
explicitly shown it is possible to renormalize the gap 
equation and effective potential in a way that is inde- 
pendent of temperature and chemical potentials. We 
have studied the phase diagram and the quasiparticle 
masses of the 0{2) x 0(2) model where the neutral kaons 
condense at sufficiently low temperature and sufficiently 
large value of the chemical potential. If the transition 
is first order, it turns out that the transition is not to 
a symmetric state but to a state with a conden- 
sate. This is in agreement with the findings of Alford, 
Braby, and Schmitt [l^. Generally, our predictions for 
the critical temperature and other quantities are differ- 
ent from theirs since our couplings do not depend on 
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the chemical potentials. Secondly, we have not made 
a high-temperature approximation of the thermal inte- 
grals in the gap equations that determine the dressed 
masses. Finally, our gap equations and effective poten- 
tial include renormalization effects. 

One drawback of the 2P1 Hartree approximation is 
that it does not obey Goldstone's theorem. We have 
shown this explicitly and quantified the deviation by the 
right-hand-side of Eq. ([70]) . For practical purposes, the 
violation is negligible which is reassuring. We therefore 
believe that the 2PI Hartree approximation is a useful 
nonperturbative approximation for systems in thermal 
equilibrium. 

The Hartree approximation and the large- iV limit are 
both mean- field approximations. It would be desire- 
able to go beyond mean field for example by includ- 
ing next-to- leading corrections in the 1/iV-expansion. In 
the context of pions, this has been done using the IPI 
1 /TV-ex pan sion , but only for vanishing chemical poten- 
tial (5ll.l52j| . Another method is the functional renormal- 
ization group method, see Refs. [H, [H, [H, [s^ and ref- 
erences therein. The functional renormalization group 
method is a nonperturbative approach that has been 
very succesful in finite-temperature field theory. The 
essence of this approach is a flow equation for the av- 
erage effective action. This flow equation cannot be 
solved exactly, but one must resort to approximations. 
A systematic approximation is the derivative expansion 
of the effective action and even the simplest truncation, 
namely the local-potential approximation, often yields 
good results. In the context of Bose condensation, func- 
tional renormalization group methods have been applied 
to the nonrelativistic case in Refs. [13, [H, [H, [13] ■ Work 
on pion and kaon condensation using these methods is 
in progress [6lj. 

APPENDIX A: RENORMALIZATION 

In this section, we first discuss the renormalization of 
the gap equations and effective potential for the 0{2N)- 
symmetric models in detail. We then briefly skecth the 
renormalization of the gap equations of the 0{2) x 0{2)- 
symmetric model. 

1. 0(2A'^)-symmetric models 

The effective potential and the gap equations contain 
divergent sum-integrals. For the theory to be renormal- 
izable, the divergent terms must be independent of the 
the temperature T and the chemical potential fi. Below 
we show explicitly that there are individual /z-dependent 
contributions to the gap equations which are divergent, 
but they cancel amongst themselves and so the coun- 



terterms are the same as those needed to renormalize 
the theory in the vacuum. 

There is a complication regarding the renormalization 
of the quartic coupling in the 2PI Hartree approxima- 
tion. The truncatation of the 2PI effective action allows 
one to define two independent four-point functions and 
thus two i ndep endent sets of counterterms associated 
with them [43|. They correspond to the two 0{N) in- 
variants SijSki and SikSji + SikSjk from which Fijki is 
built. More generally, each 0(A^)-invariant term in the 
2PI effective action has its own independent countert- 
erm. For example, the counterterms arising from the 
classical potential in Eq. (fT8|) are written as 

Isml + g^^^ (Al) 

while the counterterms arising from the term TyDq^D 
are written as 

5ml = Sm^ + ' ^ (bo , (A2) 

Smj = W + . (A3) 

The two counterterms and 5X^ and SX^ . correspond to 
the two 0(iV)-invariant terms Tr[0^]Tr[i:»] and Tr[(/)gD] 
in the expression for Tr Dq^D that scale as N and one, 
respectively. It was shown in Ref. :42:] that Suiq = Sm^ 
and(5Ao = (5A^-F2A-^. 

Renormalization of the gap equations can be done by 
applying an iterative procedure which is discussed in 
Refs. [4l|, To this end, we write the mass coun- 

terterm and the coupling constant counterterms Sm'^, 
6X^, and SX^ as power series in A: 

oo 

(5m^ = ^ Sml , (A4) 

n=l 

oo 

6X^ = (A5) 

oo 

SX^ = ^<5A^ (A6) 

n=l 

where (5rn„ is a counterterm of order A", and 6X^ and 
SX^ are counterterms of order A""'"^, respectively. Sim- 
ilarly, we write the self-energies H^ as power series in 
A: 

H, = Hf ^ + nf ) + • • • , (A7) 

where the superscript indicates the power of A. Inserting 
Eqs. d^-dMI into we can determine u'k \ Sml, 
SX^, and 6X^ by iteration. 
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The Schwinger-Dyson equation (P5|) for the three diagonal components of the propagator is 



D22 = (^0^)22+2^ [^^11+ ^22 + (2^ -2)^ i?33 
- (^o-)33 + 4 [^^ + ^22 + 2.r^^ ..33] 



(A8) 
(A9) 
(AlO) 



Since the off-diagonal parts of the self-energy vanish we have (-Do)y^ = (^)ij^ '^'-"^ * 7^ "^^"^ ^'^^^ iteration of the 
gap equation is found by ignoring the self-energy 11 on the right-hand side of the gap equations (|A8p - (|A10p . i.e. one 
replaces the full propagator by the tree-level propagator. If we denote the components of the free propagator by 
(£'o)iji we can write 



Mf 


— nil 


= (5m^ -t- 




+ 2^Af ^ 
2N ^ 


Ml 




— Srn^ + 


2N 


^2 A 
^^^2N 


Mi 




= Sm"^ + 


2N 


^2 A 
^^+2N 



)33 



(^o)ii + srf (7^0)22 + {2N - 2)jf (Do) 

JQ JQ JQ 

JQ 



T (Z?o)22 + 2iV] 

Tq JQ 



'Q 

)33 



(All) 
(A12) 
(A13) 



where we have used Eqs. (|A2|) - ()A3p . The diagram- 
matic interpretation of the procedure is shown in 
Fig. [3) The propagator in the loops are free propaga- 
tors and the first Feynman diagram on the right-hand 
side corresponds to the sum-integrals in Eqs. (jAlip - 
(jA13p . We now consider in detail the renormaliza- 
tion of Eq. (jAlip . The sum- integrals can be split into 
divergent temperature-independent and temperature- 
dependent convergent parts and so we write 

jf {Do)nn = It + II . (A14) 

I 

We begin by summing over the Matsubara frequences. This yields 



The free propagator Dq is obtained from Eq. (pS)) by 
replacing the medium-dependent masses Mi by m^. The 
first sum-integral in Eq. (jAlip then becomes 



0)11 — 



(A15) 



i 









4 








y^4/x2 -I- i(TO2 ^ ml)] + jimf — ml) 
\{m\ - ml) + 2/^2 



y^4^2 [^2+1(^2^^2)] + 1(^2 _ ^2)2^ 



COthf^) 

\2TJ 



(A16) 



We are interested in isolating the divergent parts of the sum-integral. Dropping the convergent temperature- 
dependent parts of coth (-^^) and rearranging, we obtain 



It = 



1 1 



1 1 



\{m\ ~ ml) + 2fi'^ 



^4^1^ [g2 + i(m2 + m2)] + ^{mj ~ ml)^ 



For convenience, we introduce some shorthand notation: 

El = 



(A17) 



(A18) 



14 



A 
B 



g2 + -{ml+niD 



The various factors in Eq. (jA17[) can then be compactly written as 

^1,2(9) 

i(m^ - ml) + 




-A + 



The next step is to expand wi.2('z) in inverse powers of E. The integral can then be written 



(A19) 
(A20) 

(A21) 
(A22) 



as 



1 



A~B ?,{A-B) 



2El 



8Ef 



A + B 3{A + B) 



2El 



8Ef 



2 

e'i 



A-B iiA-B) 
'^Ef^ 



8Ef 



A + B 3{A + B) 



2El 



8Ef 



-A + 3fi^ 
B 



2EI 



(A23) 



r 



We note that the integrals of the form E^ " are di- 
vergent in the ultraviolet for n < 3 and convergent for 
n > 3. Since the terms with n — 3 cancel in the integral, 
there is only one divergent term. We denote this term 
by / and it reads 



Having isolated the UV divergence, the temperature- 
independenet part of the sum-integral in Eq. (jA16|) can 
now be written as 



7^1 



1 



2 Jq ^q^ + m\ 



(A24) 



Td 



1 + 



(A26) 



This integral is easily calculated with dimensional regu- 
larization 



(47r)2 



A 



2 \ <^ 



1 



+ l + 0(e) 



(A25) where If is the finite part of /f and is given by 





■ 1 


1 

1 — 


+ 


■ 1 


1 " 






UJ2 




.^1 


W2_ 



^(mf - m^) -t- 2^2 



y^V [q2 + 1(^2 _^ „,2)j ^ 1 („,2 _ „,2)2 ^q^ + vn\ J (47r) 



where L„ = log -t- . The temperature-dependent part of the sum-integral (jA16P is given by 



h (^1 + 1) 



(A27) 



li = 



1 



1 



Wi(e'^i/^-l) tj2(e"=/'^ - 1) 
1 1 



Wi(e'^i/^-l) cj2(e'^^/'^- 1) 



;(mf — ml) + 2^2 



y^4^2 [^2 ^ 1 (^2 j„2)] _^ 1(,„2 _ ^2)2 



(A28) 
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The integrals in Eqs. (jA27|) ~ (|A28|) are now finite in three dimensions and we set e = 0. The contributions to the 
right-hand side of Eq. (jAlip from (1)0)22 and (I?o)33 are calculated in the same manner. The gap equation can then 
be written as 

Ml -m\= 5m' + ^^M±^</>2 + A| [3™? + mj + (2iV - 2)m^] 1 + 3/f + /^^ + {2N 2)11 

+3/f + + {2N - 2)/3^| . (A29) 

I 



where /| and /J are obtained by exchanging m\ for rn^ 
and vice versa in Eqs. (U27| and dMSJ, and 

^3 = -7£!l(^3 + l), (A30) 



The procedure is carried out iteratively to all orders in 
ml and A. The final result is given in Eqs. (|30l) - (|32l) 
The counterterms 5m1 5X^, and (5A^ are then expressed 
in terms of dml_i, 6X^_i, and X^_i and obtain the 
recursion relations 



^3 



1 

2^ 



dqq'^ 



+ ™i eV^^/^ 



1 



.(A31) 



The divergences in Eq. (|A29[) are now cancelled by 
choosing the counterterms appropriately. This gives 



5m\ 
5Xf 



X 



(47r)2e 

(47r)2e 

A^ 1 

(47r)2eiV 



1 + — 



1 + — 



5ml = 
^A^ = 



A 



(47r)2e 
A 

(4^ 
A 



1 

TV 



<5A:; 



5X1 



N 



(A32) 
(A33) 
(A34) 



(47r)2e N 



(A36) 
(A37) 
(A38) 



After renormalization, the first iteration of the gap equa- 
tion for All can be written as 



The counterterms are in agreement with those found by 
Fejos, Patkos, and Szep [42] ■ We note in particular that 
they are independent of temperature and chemical po- 
tential. They simplify significantly in the large- A'^ limit, 
where 5X^ = 0. Renormalizing to all orders effectively 
means that we replace the tree-level masses on the 
right-hand side of the gap equations by the medium de- 
pendent masses M,. For each /„ in Eqs. (|M7)) . (|M8)) . 
and (|MT1) we therefore define a corresponding J„, where 
the tree level masses m„ are replaced by the dressed 
masses M„. The final result is given in Eqs. ([30|) -([32 |l . 



We next consider the gap equation (|24|) . Differentiating Q with respect to (po, we obtain unrenormalized gap 
equation 



Ml - ml + ^ [31^1 +I^ + {2N- 2)11 
3I^ + I^ + {2N-2)lJ] . 



(A35) 



J 







2 2 , 

m — + 



_X_ 

2N' 



X 

2N 



JQ JQ 



D22 + i2N~2)Jl D, 



^33 



Renormalizing the gap equation (|A39p in the same way as above, we obtain 

A „ A 







2N 



(A39) 



(A40) 



We should make sure that the counterterms that renormalize the gap equations are also sufficient to renormalize 
the effective potential We next consider separately the different terms contributing to Vl. The first term is 

TrlogD"^. Taking the trace yields 



-Trlogi^-i - irf log [ujl+Col{q)] + Uf log [u:l + Coliq)] + (TV - log [u 



(A41) 
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To carry out renormalization by the iterative procedure, we expand the sum-integrals in powers of the self-energies 
11^. To first order in the self-energies, we obtain 

\Tr\0gD-' = log K + ^2(g)] ^ 1^ [^2 ^ ^2(^)] _^ _ [^2 ^ ^2(^)] 

+ 2^^ ¥gK + -?)(c.^+c.|) + 2n2 ^^^3 9^^2+,2-(A42) 



Summing over the Matsubara frequencies, each term involving log(tJ,^j -t- Lof(q)) can be written in the form 

^ \og[ul + ujj{q)\ =Xi+Xl, 
JQ 



where 



T';, - 2T / log 



1 _ p-^^{<^)|T 



(A43) 

(A44) 
(A45) 



We see right away that the temperature-dependent integrals are convergent. We now consider the temperature- 
independent term in Eq. (jA44p where n— 1,2. With the expected result in mind we rewrite wi +ijJ2 = (^1 + W2)/2 -I- 
-\- 102)12 and expand the terms in the first parantheses in E\ = \/ (p' -V nrf and the second in E2 = ^/q^ + rn^. 
This yields 



El E2 



ui[q)+uj2{q)=Ei+E2--{mi-mi) — - — _ - m^)^ IxT + 1?3 + 0{ll El) + 0{l/ E^) {MQ) 



16' 



El El 



The first two terms are the divergences we get when /i = 0. As for the other terms, writing E2 — El — {m\ — m|) 
and expanding the resulting expression in we find that the divergent parts do in fact cancel: 



'-{ml - ml) 



1 1 

El E2 



2\2 



16 



{mi - TO2) 



1 1 



-0{\/El)^0{\/El)=0'{\/El] 



(A47) 



r 



This shows that the divergent term in each sum-integral ily calculated and reads 
is on the form 



Kn = / • (A48) 

With dimensional regularization, the integral can be eas- 



{-: 


€ 


'1 


3" 




) 


e 


^2 



J 



To zeroth order in the self-energies 11^, we can write 



^Trlogl?"^ = - 



1 



1 1 



(A49) 



4(47r)2 



[ml + mi + 2{N - l)ml] - + -{It+lf) + Tli^^+ll) + {N ~ mi+ll) ,(A50) 



where 



[^^1,2(9) - -El, 2] 



"2(47r)2 



2(47r)2 



(A51) 
(A52) 
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We next consider the term TiD^^D. Using that U = D-^ - D^'^ , we obtain TiD^^D = Tr/ - TrHD, where / is 
the identity matrix. The first term vanishes in dimensional regularization and we are left with TvYlD. Since 11 is 
diagonal, we obtain 



2™ ^ 2HRTWT^+2HM 



m: 



//-2 



(A53) 



Since we are expanding the self-energies as 11^ = n'^-* + liY' + the leading term is obtained by replacing 11^ 

by n^^'' and Qi by oji. We then see that the leading term exactly cancels against the terms on the second line in 
Eq. (|A42[) . At higher orders this no longer the case. The term TrllD then gives rise to an additional term which 
after renormalization reads 



r(2) 



(A54) 



Finally, we consider the two-loop diagrams in $[-D]. The terms are all given by products of one- loop sum-integrals 
that we have already calculated 



JQ JQ JQ 



DnJl Dn 



D 



22 



Q 



D 



33 



(A55) 



Expressing this in terms of the integrals if and , we obtain 



SN 



-2 - 



(47r)2e ^ ^ 



+ It 



(4^)2e 



+ /| + /J + 2{N - 1) 



(47r)2e 



2 - 



(47r)2e 



+ 4(Ar- 1) 



(47r)2e 
(4^ 



. (A56) 



The counterterm arising from ^TrD^^D is denoted by ^TiSD^^D. Since the matrix 6Dq is diagonal, whose 

(A57) 



elements are given by Srnf, the matrix SD^ D is also diagonal. We then find 

^TtSD^'D = bmljf Dn + \5mfd D22 + {N - l)Smlj{ D. 



Again the sum-integrals are known and the result is 



(47r)2e 



(47r)2e 



/.J 



+ {N - l)6ml 



(47r)2e 



(A58) 



Adding Eqs. (|A50p . (jA56p . and (jASSp and using the counterterms in Eqs (jA32p - (IA34p . all the divergences cancel 
and we obtain the renormalized effective potential to first order 



(A59) 



The full nonperturbative effective potential is now obtained by renormalizing to all orders in the coupling constant. 
The final result is then given by Eq. (P7|) . where the integrals J^f and J'.^ are obtained from Tf and if by the 
replacements Mi. 
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2. 0(2) X 0(2)-symmetric models 



The gap equations for the dressed masses and the background field 0o foUow from Eqs. (p5| and ([6^ 



-2 [3AoZ?ll + \0D22 + AHi?33 + AhZ?44] 



= + -^4>l + -^^ [AhI?!! + \hD22 + 3A+i?33 + KDii] , 

Ml -ml = Sml + ^cbl + -^ [XnOn + XhD22 + X+D33 + 3X+Dm] , 



Mi - mi 



^"-m"+3Ao3 Dn + XnU D22 + XhJL D33 + XhJI D^^ 

Q JQ JQ 



(A60) 
(A61) 
(A62) 
(A63) 

(A64) 



Note that we have two counterterms Xq and A^ for Ao that correspond to the two 0(2)-invariant terms in the 
effective action. 

Again renormahzation can be carried out by the iterative procedure of Refs. [ill, The first step consists of 
neglecting the self-energies in the gap equations. The counterterms necessary to render them finite are given by 



1 



(1) (47r)2e 
1 



(47r)2e 



(5Ao(i) - 



1 



(47r)2e 
(47r)2e 



(47r)2e 



[2X+ml + XHml] 
{XhXo + XhX+) , 



(A65) 
(A66) 
(A67) 
(A68) 
(A69) 



where the index (1) indicates that this is the first iteration. The iterative procedure leads to the following recursion 
relations for the counterterms 



1 



(47r)2e 
1 



'^"^+(") ^ (4^ 



('^•^o(«-i) + '5^o(n-i)) rnl + <5Aff(„_i)r74 
[25A+(„_i)TO^ + 6XH(n-i)ml] , 



Ao (^2(5Ao(„_;^) + (5A^(„_i)) + Ah(5Ah(,i-i) 



SX' 



SX 



+(") 



SX 



H{n) 



(47r)2e 
1 

(47r)2e 



2A+(5A+(„_i) + -XHSXfi{n~i) 



(A70) 
(A71) 
(A72) 
(A73) 
(A74) 
(A75) 



r 



The dressed masses and M4 are equal. This follows izing the resulting expression order by order and using 
from subtracting Eq (|A63p from Eq (jA62[) and renormal- the fact that the tree-level masses m^ and m4 are equal. 
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Clearly, this holds only as long as there is no charged 
condensate. 

These recursion relations are in agreement with those 
found by Fejos, Patkos, and Szep [42] for the 0{N) x 
0(Af)-symmetry after having set N = M = 2. The 
final result for the gap equations is given in Eqs. ((65)) - 
([68]) . Finally, the counterterms needed to render the 
gap equations finite are also those needed to renormalize 
the effective potential. The final result for the effective 



potential is given by Eq. ([7T|) . 
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